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ABSTRACT 

We perform numerical simulations to study the secular orbital evolution and 
dynamical structure in the HD 69830 planetary system with the best-fit orbital 
solutions by Lovis and coworkers (2006). In the simulations, we show that the 
triplet Neptunian system can be stable at least for 2 Gyr and the stability would 
not be greatly influenced even if we vary the planetary masses from Neptune-mass 
to Jupiter-mass. In addition, we employ the Laplace-Lagrange secular theory to 
investigate the long-term behaviors of the system, and the outcomes demonstrate 
that this theory can well describe and predict the secular orbital evolution for 
three Neptune-mass planets, where the secular periods and amplitudes in the ec- 
centricities are in good agreement with those of the direct numerical integrations. 
We first reveal that the secular periods of the eccentricity ei and 62 are identical 
about 8,300 yr, i.e.. Pi = P2 — 2-k /\gi — ((71,2 are respectively, the eigenfre- 
quencies of the system), while the secular variation of 63 for the outermost planet 
has a period of ~ 18, 200 yr. Moreover, we extensively explore the planetary con- 
figuration of three Neptune-mass companions with one massive terrestrial planet 
residing in 0.07 AU < a < 1.20 AU, to examine the asteroid structure in this 
system. We underline that there are stable zones at least 10^ yr for low-mass 
terrestrial planets locating between 0.3 and 0.5 AU, and 0.8 and 1.2 AU with final 
eccentricities of e < 0.20. Still, we also find that the secular resonance vi and 
V2 arising from two inner planets can excite the eccentricities of the terrestrial 
bodies, and the accumulation or depletion of the asteroid belt are also shaped 
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by orbital resonances of the outer planets, for example, the asteroidal gaps at 
2:1 and 3:2 mean motion resonances (MMRs) with Planet C, and 5:2 and 1:2 
MMRs with Planet D. In a dynamical sense, the proper candidate regions for the 
existence of the potential terrestrial planets or Habitable Zones (HZs) are 0.35 
AU < a < 0.50 AU, and 0.80 AU < a < 1.00 AU for relatively low eccentricities, 
which makes sense to have the possible asteroidal structure in the system. 

Subject headings: celestial mechanics-methods:n-body simulations-planetary systems- 
stars:individual (HD 69830, 47 Uma, 55 Cancri, HD 160691, GJ 876) 

1. Introduction 

To date, more than 200 extrasolar planets have been discovered about the nearby stars 
within 200 pc (Butler et al. 2006; The Extrasolar Planets Encyclopaedia) mostly by the 
measurements of Doppler surveys^ and transiting techniques. The increasing numbers of the 
extrasolar planets are greatly attributed to the increasing of precision of measurement. At 
present, the observational precision has been achieved to ~ 1 ms~^ to 3 ms~^, whereas most 
observations currently have lower precision. The improvement of the observations will indeed 
induce the substantial discovery: (l)plentiful multiple systems, which ~ 20 multiple system 
involved in orbital resonance and secular interactions now have been detected; (2) much 
more low-mass companions around main-sequence stars (from Neptune-mass to Earth-mass, 
so-called super-Earths) , e.g., 55 Cancri (McArthur et al. 2004), GJ 876 (Rivera et al. 2005), 
HD 160691 (Santos et al. 2004; Gozdziewski et al. 2006);(3)a true Solar System analog, 
bearing several terrestrial planets, the asteroidal structure and the dynamical environment 
of potentially terrestrial planets in the Habitable Zones (HZs) (Kasting et al. 1993) for the 
development of life; (4)the census of the systems, and the diversity of the planetary systems 
may provide abundant clues for theorists to more accurately model the planetary formation 
(Ida & Lin 2004; Boss 2006). 

Lovis et al. (2006) (hereafter Paper I) recently reported the discovery of an interest- 
ing system of three Neptune-mass planets orbiting about HD 69830 through high precision 
measurements with the HARPS spectrograph at La Silla, Chile. The nearby star HD 69830 
is of spectral type KOV with an estimated mass of 0.86 ± O.O3M0 and a total luminosity of 
0.60 ± 0.03-^0 (Paper I), about 12.6 pc away from the Sun. In addition, Beichman et al. 
(2005) announced the detection of a large excess owing to hot grains of crystalline silicates 
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orbiting the star HD 69830 and inferred that there could be probably a massive asteroid 
within 1 AU. Subsequently, Alibert et al. (2006) and Paper I performed lots of calculations 
to simulate the system and revealed that the innermost planet may possess a rocky core sur- 
rounded by a tiny gaseous envelope and probably formed inside the ice line in the beginning, 
whereas the two outer companions formed outside the ice line from a rocky embryo and then 
accreted the water and gas onto the envelope in the sequent formation process. Hence, it 
is important for one to understand the dynamical structure in the final assemblage of the 
planetary system (Asghari et al. 2004; Ji et al. 2005; Jiang & Yeh 2006), and to investigate 
suitable HZs for hfe-bearing terrestrial planets (David et al. 2003; Fatuzzo et al. 2006; Jones 
et al. 2005; Haghighipour 2006) advancing the space missions (such as CoRot, Kepler and 
TPF) aiming at detecting them, thus one of our goals is to focus on this issue (§2.3). 

On the other hand, from a dynamical viewpoint, one may be concerned about the long- 
term stability of the system. In the present work, we further perform a large number of 
simulations of this system to examine the stability (§2.1) and to understand the secular 
evolution by means of the secular theory (§2.2). Finally, we summary the main results and 
present a concise discussion in §3. 

2. Dynamical Analysis 

In this paper, we adopt the orbital parameters of the HD 69830 system provided by 
Paper 1. In the simulations, we always take the stellar mass and the minimum planetary 
masses from Table 1, except where noted. To be specific, the mass of the host star is 
O.SGMq, and the masses of three Neptunian planets are respectively, IO.2M0, II.8M0 and 
I8.IM0, where sini = 1. We utilize N-body codes (Ji, Li & Liu 2002) to perform numerical 
simulations by using RKF7(8) and symplectic integrators (Wisdom & Holman 1991) for this 
system. In the numerical study, the adopted time stepsize is usually ~ 2% - 5% of the orbital 
period of the innermost planet. In addition, the numerical errors were effectively controlled 
over the integration timescale, and the total energy is generally conserved to 10~^ for the 
integrations. The typical timescale of simulations of the three-planet system is from 100 Myr 
to 2000 Myr. 

2.1. The Stability of the HD 69830 Planetary System 

To examine the secular stability of this system, we numerically integrated the three- 
planet system on a timescale of 2 Gyr in terms of the initial values hsted in Table 1. In 
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Figure 1, a snapshot of the secular orbital evolution of three planets is illustrated, where 
Qi = aj(l + Cj), Qi = ai{l — Cj) (the subscipt i = 1,2,3, individually, denoting Planet B, 
C and D) are, respectively, the apoapsis and periapsis distances. In the secular dynamics, 
the semi-major axis ai remains unchanged to be 0.0785 AU for 2 Gyr, whilst a2 and 
slightly undergo librations about 0.186 and 0.630 AU with small amplitudes over the same 
timescale. The variations of eccentricities during long-term evolution are followed, where 
0.05 < ei < 0.20, 0.0 < ea < 0.15 and 0.069 < eg < 0.078. In Fig.l, we note that the 
time behaviors of Qi and again show regular motions of bounded orbits for three planets, 
indicating their orbits are well separated in the secular evolution due to smaller mutual 
interactions. In the numerical study, we find the system can be dynamically stable and last 
at least for 2 Gyr, and our simulations do extend and confirm the numerical exploration of 
1 Gyr by Paper I. 

To explore the stabihty of HD 69830 with respect to the variations of the planetary 
masses, we first fix sini in increment of 0.1 from 0.1 to 0.9. In the additional numerical 
experiments, we simply alter the masses but keep all orbital parameters (Table 1), again 
restart new runs of integration of three-planet system for 100 - 1000 Myr with the rescaled 
masses. As a result, we find the system could remain definitely stable for the above investi- 
gated timescale with slight vibrations in semi-major axes and eccentricities for all planets. 
However, Paper I also pointed out that the system can even survive for the planets bearing 
the Jupiter-masses. To verify this, we further let i — 1.0°, 0.8°, 0.5°, then the planetary 
masses are in the Jupiter-mass range, and the outcomes for new integrations of 100 Myr 
show that there is no chaotic behavior occurring in this system. This may also imply that 
the stabihty of the HD 69830 system is not so sensitive to the planetary masses. 

Next, we come to the qualitative stability criteria for this system. Gladman (1993) 
analytically attained the minimum separation A between two planets of masses mi and m2 
moving on initially circular orbits required for stability. 



With the orbital data in Table 1, we then approximately estimate A12 — 0.0135 AU, for 
the couple of Planet B and C; similarly, A23 ~ 0.0460 AU, for the pair of Planet C and 
D. Note this criteria are good guidelines for stabihty but are not definitive constraints for 
nearly-circular orbits. However, Chambers et al. (1996) presented a more general stability 
criterion that shows the initial orbital separation should be no less than about 10 Rh for 



A ~ 2V3R 



(1) 



where the mutual Hill radius Rh is defined by 




(2) 
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multi-planet systems. For the Neptune-mass planets, the actual orbital separations are ~ 10 
times of the analytical values that the stability requires, while the mutual distances between 
two orbits can amount to be ~ 10 Rh-, even though these planets possess Jupiter masses. 
Therefore, from both the numerical and analytical views, we can safely conclude that the 
HD 69830 system is dynamically stable for the lifetime of the star. 

2.2. Seculeir Dynamics 

Considering a system of massive bodies of mass rrii (i — 1,N) moving about the 
central star of mass under their mutual gravitational forces and the attraction of the 
central body, if ^ m^, and these planets are with the orbits of small eccentricities and 
low inclinations, then the long-term evolution of the eccentricities and inclinations of all 
planets can be analytically described by classical Laplace-Lagrange (hereafter L-L) secular 
theory, provided that the mean motion resonance amongst them is absent. Many authors 
have adopted the secular theory to study the orbital evolution of the planets in the systems, 
e.g., Ji et al. (2003) applied this theory to explore the dynamical evolution of three giant 
planets in 55 Cancri; more recently, Adams & Laughlin (2006; and references therein) further 
investigated the dynamical interactions in the multiple planetary system by using secular 
theory and showed that the orbital eccentricities in these systems alter over secular time 
scales, of typical variations over 10^ — 10^ years. 

As for the HD 69830 system, there are three Neptune-mass planets orbiting the host 
star, assuming all planets are coplanar. Furthermore, in Table 1 we notice that the orbital 
eccentricities are smaller than about 0.10, hence one may realize that the L-L secular the- 
ory can be applied to this system. As usual, the L-L secular theory tells us the story of 
the precession of the longitudes of periastron and that of the longitudes of node, and the 
variations of the amplitudes of the eccentricities and inclinations. Herein, we simply address 
the issue of the secular evolution of the eccentricities. Bearing this in mind, the L-L theory 
is an approximation to the actual motions, because it is valid only to the second order in 
the eccentricities and inclinations and the first order in the planetary masses, however, it 
may reveal substantial clues for long-term evolution of the system. Subsequently, we will use 
the analytical results from L-L theory to compare with the numerical outcomes from direct 
integrations. 

The hnear L-L secular solutions for the eccentricities (the eccentricity eigenvectors, see 
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Murray & Dermott 1999, MD99) can be written : 

3 

Cj sin Wj = ^2 eji sin{git + (3) 

i=l 

3 

Cj COS zuj = Cji cos{git + (4) 

i=l 

where eji, Qi and /5i (j = 1,2,3, respectively for Companions B, C and D; i — 1,2,3) are 
respectively, the magnitudes of eccentricity eigenvectors, the eigenfrequencies and the phases, 
which can be calculated from the initial orbital solutions given in Table 1. In Table 2 are 
listed the magnitudes of the eccentricity eigenvectors tji evaluated from the secular theory 
by considering all three planets. Firstly, we obtain three eigenfrequencies gi = 236 ".02 yr~^, 
g2 — 79 ".76 yr^^ and g^ = 8 .49 yr^^ (see Table 2) for the HD 69830 system, and these values 
are related to the precession of the periastron of three Neptune-mass planets with periods 
of 5,491 yr, 16,248 yr and 152,650 yr, which are essentially in agreement with the results of 
Paper I. In addition, the phases are /?i = 12.12°, p2 = 244.13° and /?3 = 224.56°, and now 
all constants in the right-side of Eqs. (3) and (4) are well determined, hence we can derive 
the analytical expressions of the eccentricities ej for Planet j at any time t, and its squares 
are given by 

3 3 

= ( sm{git + (3i)f + eji cos{git + A))^ (5) 

i=l i=l 

Secondly, we estimate the period of ei in the secular orbital evolution using the L-L 
theory, from Table 2, we note that the magnitudes of the eccentricity vectors en, ei2 ^ 613, 
and it is not difficult for one to understand that the mutual interactions of two inner planets 
are much greater than that of the distant outermost planet. Hence, we have 

el ~ e^i + + 2enei2 cos{{gi - g2)t + - /^s))- (6) 

Then, the secular period of ei given by analytical method is Pi ~ 2T:/{gi — g2), about 
8,300 yr. Similarly, the secular period of 62 can be also derived, where P2 — 27r/|(g'2 ~ 
apparently, is equal to Pi. Now, we compare the analytical predictions with the numerical 
integrations. In Fig. 2, a snapshot is illustrated to perform the comparisons, while it can 
reflect the situation of total timcscalc of 2 Gyr. In this flgurc, the dash dot line denotes 
the variational eccentricities of numerical results, while the solid line represents those from 
secular theory. For the upper and middle panels. Pi and P2 are clearly evident and also 
in accordance with the corresponding long-term periods for eccentricities in the numerical 
results. In the meantime, the librating amplitudes for the eccentricities ei, 62 from the L-L 
theory {solid line), where {eimim ^imax) — (1 611 — 612], |eii-|-ei2|) = (0.063, 0.193), then ei is in 
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the range [0.05, 0.20] and in [0.0, 0.15], are well coupled with those of the numerical results 
(dash dot line) by long-term integrations. Moreover, we also notice that ei experiences to 
be maximum while 62 goes down minimum in the orbital evolution (see Fig.2), and vice 
versa, still this is because of the conversation of total angular momentum for the system. 
As aforementioned, the mutual perturbations between two inner planets are much more 
important than their interactions with the outer planet (see Table 2). On the other hand, 
in the secular evolution over 2 Gyr, oi and 02 almost keep unchanged (Fig. 1), by ignoring 
Planet D, thus we have the total angular momentum ci -^/l — + C2 a/I — e| =const (ci_2 are 
the combination constants of the planetary masses and semi- major axes), hence we obtain 
the maximum- minimum correlation of Ci and 62, de facto, such orbital evolution may again 
imply two coupled secular periods of the eccentricity in this system. 

As to Planet D, the contributions of the innermost planet (HD 69830 b) is much less 
than that of the intermediate companion (HD 69830 c), indicating the mutual interaction 
from Planet C and D is dominant, then we neglect the relevant eccentricity items arising 
from Planet B in (5), and the secular period of the eccentricity of the outermost planet is 
achieved, where P3 ~ 2n/{g2 — gs) about 18,200 yr. For the bottom panel in Fig. 2, one 
can observe the long-term eccentricity periods of P3 of the secular theory and numerical 
simulations coincide with each other. Moreover, the secular theory presents the scope of 63 
in [0.069, 0.078], which also confirms the direct numerical simulations. Lastly, we reach the 
conclusion that the outcomes from L-L theory^ are in accordance with the numerical results. 

2.3. Dynamical Architecture and Potential HZs 

Modern observations by Spitzer and HST clearly exhibit that the circumstellar debris 
disks (e.g., AU Mic and f3 Pic) are quite common in the early planetary formation (for a recent 
review see Werner et al. 2006), and Beichman et al. (2006) show that there are 13 ± 3% for 
Kuiper Belt analogs around mature main sequence stars by Spitzer programs, further they 
point out that the existence of debris disks is extreme important to the resulting detection 
of individual planets and related to the formation and evolution of planetary systems. As 
mentioned previously, Beichman et al. (2005) also provide a clear evidence of the presence 
of the disk in HD 69830, and subsequently Paper I produce the best-fit orbital solutions of 



^Adams & Laughlin (2006) studied the secular theory (MD99) by adding leading-order correction for the 
general relativity (GR) and showed that GR may act to either stabilize or destabilize the system depending 
on the specific architecture of the planetary system. To examine whether GR may play a role in the stability 
of this system, we again numerically perform new runs by modifying the equations of motion including the 
GR corrections. The outcomes show that the system still retain stable in the presence of GR. 
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three Neptune-mass planets with well-separated nearly-circular orbits, and meanwhile the 
present study confirms the secular orbital stability in this system, which may imply that 
the HD 69830 system is likely to be an analog of the Solar System, then to possess other 
terrestrial planets, and to be made of a potential asteroidal belt structure or Kuiper belt 
analog in the system. Hence, it deserves to make a detailed investigation from a numerical 
perspective. 

To further investigate the dynamical structure and potential HZs in this system, we 
performed additional simulations with the planetary configuration of coplanar orbits of three 
Ncptune-mass companions with one massive Earth-like planet. In this series of runs, the mass 
of the assumed terrestrial planet ranges from 0.01 to 1.0 M®. And the initial orbital 
parameters are as follows: the numerical investigations were carried out in [a, e] parameter 
space by direct integrations, and for a uniform grid of 0.01 AU in semi-major axis (0.07 AU 
<a< 1.20 AU) and 0.01 in eccentricity (0.0 < e < 0.2), the inclinations are 0° < 7 < 5°, 
and the remainder angles are randomly distributed between 0° and 360° for each orbit, then 
the low-mass bodies were numerically integrated with three Neptune-mass planets in the HD 
69830 system. In total, about 2400 simulations were exhaustively run for typical integration 
time spans from 10^ to 10^ yr (about 10® - 10^ times of the orbital period of the innermost 
planet). Our main results now follow. 

In Figure 3, are shown the contours of the surviving time for Earth-like planets (Fig. 
3a) and the status of their final eccentricities (Fig. 3b) for the integration over 10"'' yr, where 
the horizontal and vertical axes arc the initial a and e. Fig. 3a displays that there arc stable 
zones for the low-mass planets in the regime between 0.3 and 0.5 AU, and 0.8 and 1.2 AU 
with final eccentricities of e < 0.20. Obviously, there exist unstable zones for the nearby 
orbits of three planets, where the planetary embryos there have short dynamical surviving 
time (usually less than 10^ yr) and the eccentricities can be quickly pumped up to a high 
value ~ 0.9 (see Fig. 3b), where the evolution is not so sensitive to the initial masses. In 
fact, these adjacent bodies are in association with many of the mean motion resonances of 
the Neptunian planets and the overlapping resonance mechanism (MD99) can reveal their 
chaotic behaviors of being ejected from the system in short dynamical lifetime, furthermore 
most of their orbits are within 3Rh sphere of the Neptune-mass planets, and others are 
involved in the secular resonance with two inner companions. 

Analogous to our Solar system, if consider the middle Planet C as the counterpart as 
Jupiter, we will have the regions of mean motion resonances: 2:1 (0.117 AU), 3:2 (0.142 AU), 
3:1 (0.089 AU) and 5:2 (0.101 AU), 2:3 (0.244 AU), and in Fig. 3 we notice there indeed 
exist the apparent asteroidal gaps about or within above MMRs (e.g., 3:1 and 5:2 MMRs), 
while in the region between 0.10 AU and 0.14 AU for e < 0.10, there are stable islands 
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that the planetary embryos can last at least 10^ yr. In addition, for Planet D, most of the 
terrestrial planets in 0.50 AU < a < 0.80 AU are chaotic and their eccentricities are excited 
to moderate even high values, the characterized MMRs with respect to the accumulation or 
depletion of the asteroid belt are 3:2 (0.481 AU), 2:1 (0.397 AU), 5:2 (0.342 AU), 4:3 (0.520 
AU), 1:1 (0.630 AU), 2:3 (0.826 AU), 1:2 (1.000 AU), and our results enrich those of Paper 
I for massless bodies over two consecutive 1000-year intervals, showing a broader stable 
region beyond 0.80 AU. Note that there exist stable Trojan terrestrial bodies in a narrow 
stripe about 0.630 AU, involved in 1:1 MMR with Planet D, and they can survive at least 
10^ yr with resulting small eccentricities in the extended integrations^. The stable Trojan 
configurations may possibly take place in the extrasolar planetary systems, for example, 
Ji ct al. (2005) explored such Trojan planets moving about 47 Uma, and Gozdziewski & 
Konacki (2006) also argued that there may exist Trojan pairs configurations in the HD 
128311 and HD 82943 systems. Recently, Ford & Gaudi (2006) develop a novel method 
of detecting Trojan companions to transiting close-in extrasolar planets and argue that the 
terrestrial- mass Trojans may be detectable with present ground-based observatories, whereas 
the terrestrial Trojan planets with low eccentricity orbits about 1 AU would be probably a 
fascinating habitatable environment for other life-bearing worlds. 

Malhotra (2004) and Ji et al. (2005) showed that the secular resonances may possibly 
exist in the general planetary systems, and play a part in the dynamical distributions of 
small bodies, then shape the asteroidal architecture as well as the orbital resonances. In 
this system, if consider a planetary embryo with a mass of O.OIMq (in lunar mass), the 
related location for vi secular resonance originating from Planet B is ~ 0.321 AU, where two 
eigcnfrcquencies for the terrestrial body and innermost planet are, respectively, 232 ".33 yr~^ 
and 235". 94 yr^^, which are obtained by the L-L secular theory, as the Earth-like planet 
bears the finite mass and then it may modify the values of the in situ eigcnfrcquencies gi 
{i = 1, 2, 3) of three Neptunian planets, thus z/j ~ Qi, which also indicates that both planets 
almost share common secular apsidal precession rates in their motion. Fig. 3b exhibits that 
the eccentricity of the planetesimals could be excited by ui about 0.32 AU, and meanwhile 
the 5:2 MMR about 0.34 AU (for Planet D) drives the eccentricities to modest values. 



^It would be likely that a fraction of the stable orbits become unstable over even longer timescale, 
however these results do reveal potential asteroid structure in this system. In a statistical sense, David et al. 
(2003) investigated the stability of the Earth-like planets in the binary systems and showed that the ejection 
simulations are chaotic in the studied parameter space; in addition, Fatuzzo et al. (2006) further focus on 
the characterization of the distribution of survival times for Earth-like planets orbiting a solar-type star 
with a stellar companion and present meaningful calculations on the width of the distributions of survival 
time across parameter space, and all the studies will provide helpful information of detecting such low-mass 
terrestrial planets in future survey. 
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which may imply that the combination effects of the secular resonance and MMR can act 
as a hybrid dynamical mechanism of clearing of the small bodies in the disk. Similarly, 1/2 
resonance with respect to the intermediate companion (Planet C) occurs for a lunar-size 
planet about ~ 0.962 AU at an eigenfrequency of 76 ".72 yr~^, and the V2 region is also close 
to 1:2 resonance with the outermost Planet D, where the eccentricities for low-mass bodies 
about 1 AU can be stimulated to moderate values in the extended time integrations of 10^ 
yr. Our results are in accordance with Nagasawa et al. (2005), who pointed out that the 
secular resonances can pass through the terrestrial region from outside to inside to excite the 
eccentricities of the separated planetary embryos. Note that arises from Planet D with an 
eigenfrequency of 8". 62 yr~^, which may be at work ~ 1.488 AU. In a word, the numerical 
study may suggest that there exist asteroid structure of the HD 69830 system. 



3. Summary and Discussions 

In this work, we have studied the secular dynamics and dynamical structure of HD 
69830. We now summarize the main results as follows: 

(1) In the simulations, we show that the triplet Neptunian system can be stable at least 
for 2 Gyr and that the stability would not be greatly influenced when we shift the planetary 
masses from Neptune-mass to Jupiter-mass. The L-L theory can well describe and predict 
the secular orbital evolution for three Neptune-mass planets, where both the secular periods 
and librating amplitudes in the eccentricities are well coincidence with those obtained from 
direct integrations. We first reveal that the secular periods of the eccentricity ei and 62 are 
identical, about 8,300 yr. Account for the nature of near-circular well-spaced orbits, the HD 
69830 system may be a close analog of the Solar System. Nevertheless, as well-known in 
our Solar System, Jupiter can play a major role as a dynamical shield to bounce and throw 
away most of small bodies at times and further to prevent many asteroids from colliding 
into the Earth (or other inner terrestrial planets) over the evolution. However, if all planets 
of HD 69830 are in Neptune-mass, such protection mechanism for the potential terrestrial 
planets would be weaker and larger masses would be helpful in building a harmonic habitable 
environment. 

(2) Moreover, we extensively investigate the planetary configuration of three Neptune- 
mass companions with one massive terrestrial planet in the region of 0.07 AU < a < 1.20 AU 
to examine the existence of the potential Earth-like planets and further to study the asteroid 
structure or HZ in this system. Wc show that there are stable zones at least 10^ yr for the 
low-mass terrestrial planets locating between 0.3 and 0.5 AU, and 0.8 and 1.2 AU with final 
eccentricities of e < 0.20. Furthermore, we also find that the secular resonance vi and P2 
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arising from two inner planets can excite the eccentricities of the terrestrial bodies, and the 
accumulation or depletion of the asteroid belt arc also shaped by orbital resonances of the 
outer planets, for example, the asteroidal gaps of 2:1 and 3:2 MMRs with Planet C, and 5:2 
and 1:2 resonances with Planet D. On the other hand, the stellar luminosity of HD 69830 is 
lower than that of the Sun, thus the HZ should shift inwards compared to our Solar System. 
However, Kuchner (2003) pointed out that the habitable regime is considered to be centered 
at ttp > 1.03 AU (L/L©)^/^, such that a water planet can exist in equilibrium with stellar 
radiation, where Op is the radius of the planet's orbit and L is the luminosity of the host star, 
then Qp is approximate to 0.80 AU for HD 69830. Therefore, in a dynamical consideration, 
the proper candidate regions for the existence of the potential terrestrial planets or HZs 
are 0.35 AU < a < 0.50 AU, and 0.80 AU < a < 1.00 AU for relatively low eccentricities. 
Hopefully, the future space-based observations, e.g., CoRot, Kelper and TPF will definitely 
present a handful of samples belonging to the category of the terrestrial bodies. In final, we 
may summarize that the HD 69830 system can possess the asteroidal architecture resembling 
to the Solar System and both the mean motion resonance (MMR) and secular resonances 
(Malhotra 2004; Ji et al. 2005; Nagasawa et al. 2005) will work together to influence on 
the distribution of the small bodies in the planetary system. However, the formation of this 
system is still a mystery and needs further investigation in forthcoming study. 
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Table 1. The orbital parameters of HD 69830 planetary system 



Parameter 


Companion B 


Companion C 


Companion D 




10.2 


11.8 


18.1 


Orbital period P(days) 


8.667 


31.56 


197 


a(AU) 


0.0785 


0.186 


0.630 


Eccentricity e 


0.10 


0.13 


0.07 


Longitude of Peri. -ci7(deg) 


340 


221 


224 


Epoch of Peri.(JD) 


2453496.8 


2453469.6 


2453358 



^The parameters are taken from Paper I. The mass of the star is O.86M0. 



Tabic 2. The values of e^j (xlO^)(j = 1,2,3, respectively for Companions B, C and D; 
i — 1,2, 3) and the eigenfrequencies Qj are given from secular Laplace-Lagrange theory for 
HD 69830 planetary system, and the unit for gj is yr~^ 





1 


2 


3b 


eii 


12785986 


6481193 


387801 


e2i 


-6191399 


7506478 


677060 


esi 


67818 


-312891 


7351714 


9j 


236". 02 


79". 76 


8". 49 



*^Note that 613 and 623 terms are 
contributions from Planet D, which is 
much less than those of the mutual ef- 
fects of the two inner planets. 
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Fig. 1. — Snapshot of the secular orbital evolution of three planets is shown (for t — 10^ 
yr), where Qi = aj(l + Cj), = 0^(1 — Cj) are, respectively, the apoapsis and periapsis 
distances, and the straight lines are the semi-major axes a^. Note ai remains unchanged to 
be 0.0785 AU for 2 Gyr, and 02 and 03 slightly librate about 0.186 and 0.630 AU with small 
amplitudes. Over the long-term evolution, the eccentricity variations are 0.05 < ei < 0.20, 
0.0 < 62 < 0.15 and 0.069 < 63 < 0.078. 

Fig. 2. — Comparison of secular theory with direct numerical integration for the planetary 
system HD 69830. The eccentricity variations produced by direct numerical integration are 
shown by the dash dot lines; the corresponding eccentricity variations given by secular theory 
are shown by the solid lines. The two approaches are in good agreement and the envelope 
of the variations of eccentricity are match up extremely well. 

Fig. 3. — Left: Contour of the surviving time for Earth-like planets for the integration of 
10^ yr. Right: Status of their final eccentricities. Horizontal and vertical axes are the initial 
a and e. Stable zones for the low-mass planets in the region between 0.3 and 0.5 AU, and 
0.8 and 1.2 AU with final eccentricities of e < 0.20 (cool eccentricities). Unstable zones for 
the nearby orbits of three planets, and the planetary embryos there have short dynamical 
surviving time and the eccentricities can be quickly pumped up to a high value ~ 0.9. The 
secular resonance i/i and 1/2 can excite the eccentricities of the terrestrial bodies, and the 
accumulation or depletion of the asteroid belt are also shaped by orbital resonances of the 
outer planets, for example, the asteroidal gaps of 2:1 and 3:2 MMRs with the middle planet, 
and 5:2 and 1:2 resonances with the outmost planet. 
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